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Inelastic neutron scattering is employed to study transverse spin correlations of a Heisenberg 
S' = 1/2 chain compound in a magnetic field of 7.5 T. The target compound is the antiferromagnetic 
Heisenberg S = 1/2 chain material 2(l,4-Dioxane)-2(H20)-CuCl2, or CuDCl for short. The validity 
and the limitations of the scaling relation for the transverse dynamic structure factor are tested, 
discussed and compared to the Tomonaga-Luttinger spin liquid theory and to Bethe-ansatz results 
for the Heisenberg model. 


I. INTRODUCTION 

For interacting fermions in one dimension, the notion 
of Landau’s Fermi liquid breaks down. A new theoreti¬ 
cal framework, known as the Tomonaga-Luttinger liquid 
(TLL), takes it’s place.®® The TLL is a quantum critical 
state, since quantum fluctuations dominate down to low¬ 
est temperatures, T — > 0.® Hence, its properties follow 
universal lawsP The effective low-energy Hamiltonian is 
written in a universal form with only two key parameters, 
namely the dimensionless Luttinger parameter K , which 
quantifies the interaction strength between fermions, and 
the Fermi velocity u, which links energy and momentum 
of the low-energy dynamics. The Luttinger parameter 
value K = 1 describes free fermions, while K < 1 and 
K > 1 correspond to repulsive and attractive interac¬ 
tions, respectively.® The TLL theory fully describes the 
low-energy properties, including all thermodynamics and 
correlation functions, for a variety of entirely different mi¬ 
croscopic one-dimensional (ID) models.®^ -1 Experimen¬ 
tally, the TLL has been realized in quasi-ID conductor 
material^®, fractional quantum Hall fluids^®®, or quan¬ 
tum wires, e.g. carbon nanotubes 12 13 . 

Some of the best realizations of the TLL are found in 
certain ID quan tum spin systems, including spin chains 
and ladders!^® Since spins are the relevant degrees of 
freedom in this case, they are referred to as Tomonaga- 
Luttinger spin liquids (TLSL). The humble antiferromag¬ 
netic (AF) Heisenberg ,5 = 1/2 chain is one of the old¬ 
est and best known examples. Real materials described 
by this model have provided a convenient platform for 
the study of TLSL physics using e.g. thermodynamic 
measurements ESEH neutron spectroscopy P®®^ or nuclear 
magnetic resonance (NMR) experiment J 22 l 23 l One at¬ 
traction of spin systems is that their fermionic interac¬ 
tions are continuously tunable by applying an external 
magnetic held. This feature has been successfully ex¬ 
ploited in several recent experimental studies EHHUl For 
the simple AF Heisenberg S = 1/2 chain, the Luttinger 
parameter is expected to vary continuously from K = 0.5 
in zero held to K = 1 at saturation.® This implies that 


the low-energy spin correlation functions, which are uni¬ 
versally expressed in terms of K and it, become field- 
dependent as well. Neutron scattering experiments di¬ 
rectly probe these correlation functions, and can thus be 
used to study the evolution of the TLSL in applied helds. 

To date, neutron scattering measurements of TLSL 
correlation functions in AF Heisenberg S' = 1/2 chains 
have only been performed in the absence of a magnetic 
held pStSA] Fhe ma j n result is that the AF dynamic struc¬ 
ture factor S(ir, uj) and the local dynamic structure factor 
S(u>) can be written as universal scaling forms of w/T, 
with the scaling functions explicitly given by TLSL the¬ 
ory with I\ = 0.5. As discussed below, extending such 
measurements to the case of applied magnetic helds and 
consequently larger K values, presents unique experimen¬ 
tal challenges. In the present work we attempt to over¬ 
come these difficulties. We perform hnite-temperature 
measurements of in the spin chain compound 

2(l,4-Dioxane)-2(H 2 0)-CuCl2, or CuDCl for short, in a 
magnetic held of roughly half the saturation value. De¬ 
spite a limited dynamic range, we hnd that the mea¬ 
sured dynamic structure factor does indeed show hnite- 
temperature scaling in agreement with TLSL theory. 


II. EXPERIMENTAL CONSIDERATIONS 

In the absence of an external magnetic held, probing 
TLSL physics in S’ = 1/2 Heisenberg AF materials using 
inelastic neutron scattering is comparatively straightfor¬ 
ward. At low temperatures, the excitation spectrum is 
a multi-spinon continuum with a sharply dehned lower 
bounds However, only a part of these correlations are 
actually described by the TLSL theory,®^ since it relies on 
a linear dispersion of excitations. Figure [lja) shows the 
spinon continuum and the approximate region in energy- 
momentum space, where the TLSL predictions for the 
dynamic structure factor are applicable. This region is 
relatively large, leaving a wide dynamic range available 
for investigations. The spin structure factor is generally 
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polarization-dependent, being defined as: 

(q ’ w) = 2^ / < 5 “ (x ’ *) (°»°)> e- i(qx - wt) dtdx. 

(1) 

In zero applied field, spin fluctuations are isotropic and, 
thus, S aa (q, ui) is the same for all spin components a. 
Therefore, there is no need to discriminate between dif¬ 
ferent polarization channels in a neutron-scattering ex¬ 
periment, which greatly facilitates the measurements. 

In a nonzero external magnetic field, the constraints 
on probing TLSL physics with neutron scattering are 
much more severe. A field along the z axis breaks the 
full rotational symmetry, so that the spin correlations 
become anisotropic, S xx (q, w) = S yy (q, to) A S zz (q, w). 
The total spectrum now consists of distinct and over¬ 
lapping longitudinal and transverse contributions. As 
shown in Figure [ljb) , each of these is a continuum. The 
corresponding sharp lower bounds are now dis tingu ish- 
able, with incommensurately positioned minimaP'^J The 
incommensurability is directly proportional to the field- 
induced magnetization, so that the minima of the spec¬ 
trum are shifted from the zero-field commensurate posi¬ 
tions by Sq ~ 2t r (S z ) = % ££Li (Sf) at T -> O.^Such 
a spectacular restructuring of the excitation spectrum 
has been confirmed experimentally by neutron-scattering 
studies, including those on the well-known spin-chain ma¬ 
terial CuPzNp3 In application to the present problem, 
the consequences of this spectral changes are three-fold. 
First, the lower bound of either continuum is lowered 
compared to the one in zero field and, thus, the linear- 
dispersion region where the TLSL notion may apply is re¬ 
duced. Secondly, the TLSL theory predicts different scal¬ 
ing forms for the longitudinal and transverse components 
near their respective minima. While it is in principle pos¬ 
sible to discriminate between different polarization chan¬ 
nels using polarized neutrons, in practice the correspond¬ 
ing intensity penalty makes the experiment exceedingly 
complicated or even infeasible. The only solution is to 
further shrink the measurement window to totally avoid 
the overlap between continua of different polarizations, 
as illustrated in Figure 0b). Thirdly, only considering 
one polarization channel at least halves the net intensity 
compared to the zero-field case where the scattering of 
all polarization channels can be exploited. 

One way to overcome these obstacles is to increase the 
energy range by using a target spin-chain compound with 
a larger exchange constant J. However, substantially 
magnetizing such a material would require unattainably 
large magnetic fields. The alternative is to use a com¬ 
pound with a small J, but to collect the data with tight 
energy and momentum resolutions. Unfortunately, res¬ 
olution always comes at the expense of intensity. The 
answer to this dilemma is to make use of very large single¬ 
crystal samples, ideally as large as the neutron beam it¬ 
self. This is the approach adopted in the present study. 

Guided by the potential of growing very large single 
crystals, for our experiments we selected the Heisen¬ 
berg S = 1/2 chain material CuDClP^ It has a mono- 
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Figure 1. A sketch of the spin excitation spectrum of an AF 
Heisenberg S = 1/2 chain in zero magnetic field (a) and with 
magnetic field applied along the a axis, as in the experiment 
(b). The numerical values correspond to the values of the 
target material of the present work, CuDCl. The approximate 
region where correlations are governed by TLSL physics is 
hatched. 


clinic structure, space group C2/c, and lattice constants 
a = 17.43A, b = 7.48 A, c = 11.82 A, and /3 = 119.4° 
at 173 KpS The spin chains run along the crystallo¬ 
graphic c axis and the leading term in the Hamiltonian is 
the nearest-neighbor Heisenberg exchange with exchange 
constant J ~ 0.9 meVpS The chains are formed by Cu 2+ 
ions, which are surrounded by two 0 2 ~ and Cl~ ions 
(Fig# The resulting C1-0-C1-0 plaquettes are slightly 
tilted with respect to the neighboring ones in the same 
chain, resulting in a staggered arrangement. The indi¬ 
vidual chains are well separated from each other by 1,4- 
Dioxane molecules along the a and b axes, so that the 
interchain interactions are expected to be very small. To 
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date, no magnetic ordering has been reported in this sys¬ 
tem. 



Figure 2. Crystal structure of CuDCl. The magnetic Cu 2+ 
ions are linked by Cl - ions and water molecules and form 
chains along c with coupling constant J. Individual chains 
are separated by 1,4-Dioxane molecules. The crystallographic 
unit cell is indicated by dotted lines. The Cu 2+ ions are sur¬ 
rounded by C1-0-C1-0 plaquettes that are slightly tilted rela¬ 
tive to one another in a staggered arrangment. The strongest 
interchain coupling, J ±, is expected to be along the b axis. 


III. EXPERIMENTAL DETAILS 

All experiments were performed on fully-deuterated 
single-crystal samples of CuDCl. The crystals were 
grown in a nitrogen atmosphere from a methanol solu¬ 
tion containing water, anhydrous CuCl 2 and 1,4-Dioxane. 
The crystals are highly unstable in ambient air. The 1,4- 
Dioxane molecules evaporate and leave CuCl 2 powder if 
the CuDCl crystal is not stabilized by an overpressure of 
1,4-Dioxane in the atmosphere or by an applied pressure. 

Specific-heat data were collected on a commercial 
Quantum Design physical property measurement system 
(PPMS) using a Quantum Design dilution insert for mea¬ 
surements at lowest temperatures down to 50 mK. The 
sample was covered with Apiezon N grease in order to 
prevent the 1,4-Dioxane molecules from evaporating. Be¬ 
fore the sample was mounted, the specific heat of the used 
Apiezon N was measured in order to subtract its contri¬ 
bution from the total specific heat. Each data point was 
measured with a temperature rise of 2 % of the sample 
temperature. 

Inelastic-neutron-scattering experiments were per¬ 
formed on the time-of-flight spectrometers OSIRIS at 
ISIS and DCS at the NIST Center for Neutron Research 
(NCNR). Fully-deuterated single crystals of mass 11 g 
(on OSIRIS) and 12.5 g (on DCS), respectively, were 
used. The samples were sealed in He atmosphere and 


a small amount of deuterated 1,4-Dioxane was added to 
generate an overpressure of 1,4-Dioxane inside the sample 
containers. The crystals were aligned with the incident 
neutron beam in the ac plane. A 7.5 T vertical magnet 
(H || b) was employed on OSIRIS and a 10 T vertical 
magnet was used on DCS. Both experiments were per¬ 
formed with dilution inserts for the magnets. The final 
neutron wavelength was set to 6.66 A on OSIRIS. On 
DCS an incident neutron wavelength of 8.00 A was used. 
The experimental energy resolution, defined as the full 
width at half maximum (FWHM) of the elastic incoher¬ 
ent scattering was 0.03 meV on OSIRIS and 0.06 meV 
on DCS. 


IV. EXPERIMENTAL RESULTS 

Specific heat of CuDCl was measured with a magnetic 
field applied perpendicular to the chain direction. Re¬ 
spective data for a magnetic field along the a* and b 
axes are displayed in Fig. [3] In zero magnetic field 
three features can be clearly discerned. At high tem¬ 
peratures (T > 10 K) the contribution from phonons is 
dominant. The data below 10 K, in particular the max¬ 
imum at about 3 K, are in very good agreement with 
the theoretical results for an AF Heisenberg 5 = 1/2 
chain from Johnston et all ^ with J = 0.92 meV as the 
only free parameter (solid line in Fig. § P The constant 
low-temperature part of the C/T curve spans over more 
than one decade in temperature and is characteristic of 
the TLSL phase. Below 100 mK, at the lowest acces¬ 
sible temperatures, the specific heat increases abruptly, 
suggesting the onset of three-dimensional magnetic long- 
range order. 

In applied fields, the ordering peak moves to higher 
temperatures, which is consistent with an expected in¬ 
crease of the ordering temperature Tn in weakly cou¬ 
pled spin chains. 31 33 However, a new maximum emerges 
around 0.3 K which is not a feature of ideal ID Heisen¬ 
berg 5=1/2 antiferromagnetsP At higher fields, this 
new maximum merges with the one initially seen around 
3 K. 

The neutron spectra of CuDCl were first measured on 
both spectrometers in zero magnetic field at base tem¬ 
perature. Figure |4|a) is a false-color plot of the neu¬ 
tron intensities collected on DCS in the vicinity of the 
ID AF zone-center at 85 mK. These data were taken at 
several sample orientations, then projected along the a* 
axis onto the relevant ( l , hui) plane. The result agrees 
well with previous neutron studiesP 

The spectra in a magnetic field of 7.5 T were collected 
at 2.0 K and 4.0 K on OSIRIS, as well as at 0.085 K, 
0.25 K, 0.80 K and 1.4 K on DCS. The data taken at 
low temperatures in zero field were used to determine 
the background for the high-field experiments. For each 
energy transfer, a constant background was fitted to the 
part of the zero-field data which is not affected by the 
continuum of the magnon excitations, i.e. which lies 
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Figure 3. Symbols: measured specific-heat of CuDCl for var¬ 
ious magnetic fields applied perpendicular to the chain direc¬ 
tion. An offset of 0.5 J/mol/K 2 , illustrated by an arrow in the 
graph, is added between the data sets for visibility. The solid 
line is the theoretical result for an AF Heisenberg 5 = 1/2 
chain in zero field with J = 0.92 meV, as explained in the 
text. Dashed lines are DMRG data taken from reference EH 
for a Heisenberg 5=1/2 chain with a staggered magnetic 
field of approximately 0.1 T. 


outside of the boundaries indicated by the solid lines 
in Fig. 0a). The resulting background-subtracted spec¬ 
trum for 7.5 T and 85 mK obtained on DCS is depicted in 
Fig. 0b). As expected for Heisenberg 5 = 1/2 chains, 
the spin correlations parallel and transverse to the ap¬ 
plied magnetic field show minima at different wave vec¬ 
tors. Based on the known values for the coupling con¬ 
stant J and the magnetic field H , and assuming a g 
factor g = 2, the expected boundaries of the continua 
for transverse and longitudinal excitations are shown as 
solid and dashed lines, respectively. The longitudinal 
excitations have an energy minima at l = — 1 ± 0.27, 
which corresponds to q = n ± 0.85 for the structure of 
CuDCl. From this figure it can be seen that at the ID 
AF zone center, l = —1, at low energies, one only mea¬ 
sures the transverse spin correlations, corresponding to 
S xx (tt,u}). Integrating the experimental data in a nar¬ 
row range —1.025 < l < —0.975 [dotted rectangle in 
Fig. 0b)] yields the measured energy dependence of this 
quantity, summarized in Fig. 0for different temperatures 
and both experiments. 

From Fig. 0b) it is immediately obvious that an ap- 
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Figure 4. False-color representation of the inelastic neutron 
spectra measured in CuDCl on DCS at 85 mK in zero mag¬ 
netic field (a) and at 7.5 T (b). The solid and dashed lines 
represent the lower boundaries of the continua for transverse 
(solid) and longitudinal (dashed) excitations in an ideal AF 
Heisenberg 5 = 1/2 chain compound with J = 0.92 meV at 
7.5 T. The dotted rectangle indicates the region of interest 
for the present study. The left and the right boundary of 
the rectangle mark the limits of the integration for obtaining 
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plied magnetic field opens an energy gap in the spin ex¬ 
citation spectrum, even though at zero field it is gap¬ 
less. The energy gap is approximately A ~ 0.18 meV 
at 7.5 T. The behavior is very similar to the one pre¬ 
viously studied in the spin-chain compound CDCP^ 37 ! 
In that case, it was attributed to a staggered spin field 
due to a combination of a uniform magnetic field and a 
staggered gyromagnetic tensor in the spin chains. 16 E^SH 
For CuDCl, a staggering of the g tensor would be a nat¬ 
ural consequence of the above-mentioned staggered tilt¬ 
ing of the Cl-O-Cl-O plaquettes. A staggered g tensor 
also accounts for the observed double-maximum shape 
of the C/T curves in CuDCl. The measured data are 
in good qualitative agreement with previous density- 
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Figure 5. Dynamic structure factor at the AF wave vector, 
l = —1, versus the temperature-scaled energy hu/(k.BT) for 
all neutron-scattering data obtained on OSIRIS and DCS at 
7.5 T for hw < 0.42 meV. The energy gap is clearly visible 
below 1 K as a drop of intensity in the low-energy tail of 
S (7r, lo). The error bars correspond to the statistical errors of 
the neutron-scattering measurements. 


matrix renormalization-group (DMRG) calculations for 
the Heisenberg S = 1/2 chain in a staggered magnetic 
field (dashed lines in Fig. 0.® 


V. DISCUSSION 


The ideal Heisenberg 5 = 1/2 chain below saturation. 
being in the TLSL state, is always quantum critical.^ 
However, the opening of the gap A in the excitation spec¬ 
trum is a serious complication in the quest to study quan¬ 
tum critical properties in CuDCl. Indeed, the symmetry¬ 
breaking staggered field produced in CuDCl by a uniform 
applied field takes the system away from the TLSL crit¬ 
icality. But even in this case, TLSL physics may be ac¬ 
cessible in the quantum critical regime at elevated tem¬ 
peratures or on short time scales P* Specifically, quantum 
critical fluctuations dominate for 


( 2 ) 


With this in mind, in all the subsequent analysis we only 
used those data points for which either the temperature 
or the energy transfer are much larger than the gap A, 
i.e. 


(T > 2.8 IC) or (Hu > 0.24 meV), (3) 

where the numerical limits are chosen as the measured 
energy resolution on DCS added to the gap energy. An 
additional restriction is set by the required linearity of 
the dispersion, as discussed in section |TT] In the present 
case for 7.5 T [see plot of the lower continuum boundary 
in Fig. 0b) and Fig. 0b)], the upper limit of the energy 
transfer and the temperature was chosen to be: 

k B T, Hu < 0.4 meV. (4) 


Finally, another effect which can potentially drive the 
system away from TLSL quantum criticality is three- 
dimensional interchain coupling. Fortunately, it is neg¬ 
ligibly small in CuDCl for the temperatures and energy 
transfers defined above (k B XN| 7 . 5 T ~ 0.02 meV <C A, 
cf. Fig. 0. Indeed, from the zero-field ordering tem¬ 
perature determined from the specific-heat data, In 
55 mK, the interchain coupling can be estimated to be 
Jl < 4 ^eV<C A, following reference 1351 

In general, the correlation function at the critical wave 
vector q = 7r is expected to have the following scaling 
form in the quantum critical regime: 

5(tt,u;)=T- q $(|), (5) 


with the scaling exponent a and the scaling function $. 
Quantum criticality implies that the microscopic Hamil¬ 
tonian parameters J or H do not enter this expression 
explicitly. The only relevant energy scale is set by the 
temperature itself. While the narrow dynamic range and 
the small number of different measurement temperatures 
prevent the confirmation of such a scaling law in a model- 
free fashion - as was done in our recent studies on mag¬ 
netized spin ladders,® or anisotropic spin chains at the 
Ising quantum critical point,® where the experimental 
considerations were less restrictive - the CuDCl data are 
sufficient to verify the specific prediction for the scaling 
function <f> and the scaling exponent a provided by the 
TLSL theory®!] 
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where T is Euler’s gamma function. 

Using only those S(n,u) data points that satisfy the 
above-mentioned conditions on temperature and energy 


(k B T A) or (Hu A). 
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Figure 6. Scaling plot for the Luttinger parameter K = 0.64. 
The solid line shows the universal scaling function given in 
eq. §. The error bars correspond to the statistical errors of 
the neutron-scattering measurements. 


transfer, Eqs. (§, (§ and ([7]) were fit to the data. The 
two only adjustable parameters were an overall constant 
scale factor and the Luttinger parameter K. The best fit 
is obtained using I\ = 0.64(5). The result is shown as 
a solid line in the scaling plot in Fig. [6] It can be seen 
that the data collected at different temperatures indeed 


align on a single curve with this choice of the scaling ex¬ 
ponent. Therefore, the validity of the universal scaling 
behavior of transverse spin correlations is verified. The 
measured value of K is clearly larger than K = 0.5 for 
the unmagnetized chain. Furthermore, it is in excellent 
agreement with the Bethe-ansatz result for AF Heisen¬ 
berg 5 = 1/2 chains with J = 0.92 meV at 7.5 T, which 
predicts K ss 0.65.® 

VI. CONCLUSIONS 

Despite complications intrinsic to the selected proto¬ 
type material, we have demonstrated that high-resolution 
neutron spectroscopy provides a means of probing the 
field-evolution of critical TLSL correlations in a partially 
magnetized Heisenberg 5=1/2 chain antiferromagnet. 
Furthermore, it was shown that the Luttinger parame¬ 
ter K is increased in a magnetic field compared to the 
value obtained in the absence of a field in accordance 
with Bethe-ansatz calculations. 
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